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The loop quantum cosmology "improved dynamics" of the Bianchi type IX model 
are studied. The action of the Hamiltonian constraint operator is obtained via 
techniques developed for the Bianchi type I and type II models, no new input is 
required. It is shown that the big bang and big crunch singularities are resolved 
^^ ■ by quantum gravity effects. We also present the effective equations which provide 

. modifications to the classical equations of motion due to quantum geometry effects. 
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INTRODUCTION 



cr. 

;h ! Loop quantum cosmology (LQC) [1, 2] is an approach to quantum cosmology following 

the ideas of loop quantum gravity (LQG) [3-5]. One of the major results of LQC in the 
homogeneous and isotropic Friedmann- Robertson- Walker (FRW) models is that, while gen- 
eral relativity approximates the dynamics very well in the low (with respect to the Planck 
lO \ scale) curvature regime, the classical big bang singularity is avoided: when the matter en- 

1^ ■ ergy density approaches the Planck energy density, deviations from general relativity become 

ly-^ . significant and a "quantum bounce" due to quantum gravity effects occurs when the matter 

iX ■ energy density reaches a critical energy density of the order of the the Planck density [6-14]. 

O \ More recently, it has been shown that the singularity is also resolved in the improved dy- 

^-^ ■ namics approach of loop quantum cosmology in the anisotropic Bianchi type I and type II 

cosmological models [15, 16] and in the hybrid loop-Fock quantization of the inhomogeneous 
Gowdy model [17]. The goal of this paper is to extend the LQC improved dynamics analysis 
r> ' of the Bianchi type I and type II models to the more complicated Bianchi type IX models, 

c^ ' At the classical level, the Bianchi IX model has a much richer phenomenology than 

Bianchi I and II models as it displays Mixmaster dynamics as the singularity is approached 
[18]. In essence, a space-time which exhibits Mixmaster dynamics is one which can be de- 
scribed for long periods of time (known as epochs) as a Bianchi I space-time characterized 
by three anisotropic expansion rates. Such a space-time will occasionally undergo a "Mix- 
master bounce" from one epoch to another where the three expansion rates change in a 
specific manner. Bianchi I models approach the singularity in a rather straightforward way 
as they do not undergo any Mixmaster bounces while Bianchi II models may undergo a 
single Mixmaster bounce between two epochs as the singularity is approached (see [18] and 
references therein). The Bianchi IX model, on the other hand, undergoes many Mixmaster 
bounces and this behaviour is chaotic [18, 19]. Since much of this behaviour occurs when 
the curvature is of the Planck scale, quantum gravity effects cannot be neglected and the 
Mixmaster behaviour may be significantly modified when they are taken into account. 
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Bianchi IX models are also thought to play an important role near generic singularities 
in classical general relativity. The Belinskii, Khalatnikov, Lifshitz (BKL) conjecture sug- 
gests that as a generic space-like singularity is approached, time derivatives dominate over 
spatial derivatives whence physical fields at each point evolve independently from those at 
neighbouring points. Dynamics can therefore be approximated by the ODE's used in ho- 
mogeneous space-times, most generally a Bianchi IX solution with a massless scalar field 
[20, 21]. Since other matter fields do not contribute to the dynamics as the singularity is 
approached, we will only consider the case of a massless scalar field in this work. There has 
recently been a considerable amount of numerical work supporting this paradigm (see, e.g., 
[18]) and the conjecture has also been rewritten in terms of variables suitable to a loop quan- 
tization [22]. If the BKL conjecture is indeed correct, it follows that a good understanding 
of the quantum dynamics of the Bianchi IX model in the deep quantum regime could lead 
to major insights into the behaviour of generic space-times in regions where the curvature 
reaches the Planck scale. 

Because of the Bianchi IX model's importance, it has already been the subject of studies 
within the framework of loop quantum cosmology, both in a pre-/io-type Hamiltonian frame- 
work [23, 24] and in a spin-foam-inspired dipole cosmology model (first introduced for the 
isotropic case in [12]) which also allows inhomogeneities [25]. However, it is important to 
study the improved /ij-type dynamics of LQC since it has been shown that the predictions 
of the pre-/io approach are unphysical in the infrared limit. In particular, in isotropic cosmo- 
logical models quantum gravity effects can become important at energy densities arbitrarily 
below the Planck scale in this scheme. To ensure that quantum gravity effects only become 
important at the Planck scale, one must instead use the improved dynamics approach in the 
Hamiltonian framework^. On the other hand, since the dipole cosmology model presented 
in [25] is inspired by spin foam models, that approach is complementary to ours and it will 
be interesting to compare the results of these two frameworks. 

As pointed out above, chaotic behaviour appears as the singularity is approached in 
classical Bianchi IX space-times. It has been argued in the pre-/io LQC treatment of the 
model that this behaviour is avoided in LQC due to quantum gravity effects [24]. In essence, 
the argument is that quantum gravity effects become important before a significant number 
of Mixmaster bounces occur. Since the quantum gravity effects are repulsive, the space-time 
will exit the Planck regime having only undergone a small number of Mixmaster bounces 
and hence the dynamics are not chaotic. In this paper, we see that in some cases this occurs 
already in the effective theory which incorporates the quantum geometry effects into the 
dynamics. However, we cannot yet show that this is a generic result. We will go into more 
detail in Sec. IV. 

The outline of the paper is as follows. In Sec. II, we will briefly review the necessary 
classical properties of the Bianchi type IX model in order to proceed with the quantization. 
In Sec. HI we will study the quantum properties of the model, flrst recalling the kinematics 
which are the same as for the Bianchi type I and type II models studied in [15, 16]. We will 
then study the Hamiltonian constraint operator for the Bianchi IX model with a massless 
scalar fleld as the matter fleld. The Hamiltonian constraint operator gives an evolution 
equation where the scalar fleld acts as a relational time parameter. The dynamics of the 
model are obtained by using the same technology that was developed during the study of the 



^ It is possible that once the curvature reaches the Planck scale a scheme other than /i^ may be the correct 
one but we will only consider the fli scheme here. 



improved LQC dynamics of the Bianchi type I and type II models [15, 16]; there is no need 
to introduce any new operators. In Sec. IV we will derive effective equations which provide 
the first order quantum corrections to the classical equations of motion and in Sec. V we 
summarize our results and discuss open issues. 

II. CLASSICAL THEORY 

In Bianchi models [26-28] , one restricts oneself to those phase space variables which admit 
a 3-dimensional group of symmetries which act simply and transitively. The symmetries 
allowed in the Bianchi IX group are the three spatial rotations on S^. It follows that the 
three Killing (left invariant) vector fields ^" satisfy^ 

[L ij] = -e\,L (2.1) 

where the structure constants are given by the completely antisymmetric tensor lijk times 
2/ To where Tq is the radius of the 3-sphere with respect to the fiducial metric, lijk is defined 
such that 6123 = 1) note that the internal indices i, j, A;, . . . can always be freely raised and 
lowered. There is also a canonical triad e" — the right invariant vector fields — which is Lie 
dragged by ^". It is convenient to use e" and its dual co-triad w^ as fiducial frames and 
co-frames. They satisfy: 

2 1 • • 

[ci, ej\ = e'^ij 4, du'' = — e^ijUj' A Cj^ . (2.2) 

The form of the equations above indicates that M admits global coordinates a G 
[0, 27r), /3 G [0, tt) and 7 G [0, Atx) such that for Tq = 2 the Bianchi IX co-triads have the form 

uj\ = sin/3sin7(da)a + cos7(d/3)a, 

ujI = -sin/3cos7(da)a + sin7(d/3)a, (2.3) 

ul = cos f3 (da) a + {d-f). 

The fiducial co-triads determine a fiducial 3-metric qab '■= <^a<^bii 

Qab^x'^dx^ = da^ + d/32 + d-f^ + 2cosf3dad-f. (2.4) 

It follows that ^/q = sin /3 and one can see that qab is the metric of a 3-sphere with a volume 

Vo = 167r^, this agrees with Vo = 27i'^r^ for To = 2 as specified above. Finally, we introduce 

1/3 
the length-scale io = Vo for later convenience. 

In diagonal Bianchi models, the physical triads e" are related to the fiducial ones by^ 

oo: = a\T)u;: and a,(r)e^ = e^ (2.5) 
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Here we are following the convenctions used in [25]. A different, although equivalent, choice is used in 
[10, 23, 24] where the structure constants differ by an overall sign. 

There is no sum if repeated indices are both covariant or contravariant. As usual, the Einstein summation 
convention holds if a covariant index is contracted with a contravariant index. 



where the a^ are the three directional scale factors. 

For later use, let us calculate the spin connection determined by physical triads e". Since 
r^ is given by 

r^ = -e''^ e] (dyaUJ,^, + ]-eluj[dy,uj,^^ , (2.6) 
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it follows that 



r; = £^-^-^U-. (2,7) 

To \a2a3 03 02/ 

where e := €123 is +1 for right-handed physical triads and —1 for left-handed physical triads. 
Note that while the lijk appearing in the Bianchi IX structure constants are not affected 
by the handedness of the physical triads, eijk and e on the other hand are affected by the 
handedness of e^. F^ and F^ can be obtained by permutations of Eq. (2.7). 

As is usual in LQC, we will now use the fiducial triads and co-triads in order to introduce 
a convenient parametrization of the phase space variables E^ and A^^. Because we have 
restricted ourselves to the diagonal model and these fields are symmetric under the Bianchi 
IX group, from each equivalence class of gauge related phase space variables we can choose 
a pair of the form 

Et = ^^^qet and A^ = ^ul, (2.8) 

where, as spelled out in footnote 3, there is no sum over i. Note that the length io plays a 
similar role to that of the lengths of the fiducial cell in noncompact space-times in terms of 
the form of the basic variables {A, E). In this case the manifold is compact and there is no 
fiducial cell. 

It is straightforward to relate the scale factors a^ to the pji 

pi = sgn(ai) 10203 1^0 > P2 = sgn(a2)|aia3|^o > Ps = sgn(a3)|aia2|£o , (2.9) 



it follows that \/\q\ = \/ \piP2P3\V^^ \^. 

Thus, a point in the phase space is now coordinatized by six real numbers {pi,d). One 
can use the symplectic structure in full general relativity to induce a symplectic structure 
on the six-dimensional phase space. The non-zero Poisson brackets are given by 

{c^, pj} = 8nG^ 5ij , (2.10) 

where 7 is the Barbero-Immirzi parameter. 

Our choice (2.8) of physical triads and connections has fixed the internal gauge as well 
as the diffeomorphism freedom. Furthermore, it is easy to explicitly verify that the Gauss 
and the diffeomorphism constraints are automatically satisfied due to Eq. (2.8). Thus we 
are left with the Hamiltonian constraint 

Ch= / ^: L- e'hiFj-{l+l')nab'')+NH^^tt d^x^O, (2.11) 

Jm '-167rG7 \/k/ ^ ^ -' 



where Fab'' and Qab'^ are the curvature of A^ and F^ respectively, while "Hmatt is the matter 
Hamiltonian density. The ~ indicates that Ch is a constraint and must vanish for physical 
solutions. Since we are most interested in the gravitational sector, our matter field will 
consist only of a massless scalar field T which will later serve as a relational time variable 



a la Liebniz. (Additional matter fields can be incorporated in a straightforward manner, 
modulo possible intricacies of essential self-adjointness.) Thus, 

Since we want to use the massless scalar field as relational time, it is convenient to 
use a harmonic-time gauge, i.e., assume that the time coordinate r satisfies Or = 0. The 
corresponding lapse function is A^ = \/\piP2P3\- With this choice, the Hamiltonian constraint 
simplifies considerably. 

In terms of p,, the first component of the spin connection is given by 

To V Pl P3 P2j 

the other two spin connection components can be obtained via permutations. The curvature 
of r* is in turn 
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= I f 3^ + 2 J1 - 2^ - 2^ - ^ - ^"l ulu^, (2.14) 

ri V Pl P2P3 P3 P2 PI P2 J 

the other components of Vtab' can again be obtained via permutations. 

Finally, it is straightforward to calculate the curvature of A\. For example. 

Using these results, one finds that the Hamiltonian constraint (2.11) is given by 

1 ( 2eoe. . 

^H ^ ~ Q^r<^.2 I P1P2C1C2 + P2P3C2C3 + P3P1C3C1 H [P1P2C3 + P2P3C1 + P3P1C2) 




2A + 2pl + 2pl~mi)'(?^)'(?^ 



(2.16) 



Note that the constraint for the closed isotropic case is recovered for pi = P2 = P3 while 
the Bianchi I constraint is recovered in the limit Tq — )■ 00 or, equivalently, do — ?■ 0. We will 
take advantage of this correspondence and set ro = 2 for the remainder of the paper. The 
Bianchi I limit can be obtained by taking £0 — ^ 0. 

One can now derive the time evolution of any classical observable O by taking its Poisson 
bracket with Ch'- 

d = {0,Cu}^ (2.17) 

where the 'dot' stands for derivative with respect to the harmonic time r. This gives 

Pi = — [ P2C2 + P3C3 + Le , (2.18) 

7 V Pl J 



Cl = --(p2ClC2+P3ClC3 + £.eb2C3+P3C2)+^^(l + f)fpi + ||-^-^H. (2.19) 

7 \ V 2pf 2pi 2pi J J 

As usual, the other equations of motion can be obtained by permutations. Any initial data 
satisfying the Hamiltonian constraint can be evolved by these equations of motion. It is 
particularly interesting to study the Hubble rates Hi which are given by 

H. = ^^. (2.20) 

tti at 

where t is the proper time and is related to the harmonic time r (which is the time coordinate 
used until now) by 

' ' ' (2.21) 



It follows that the Hubble rates are related to the {ci,pi) by, e.g.. 



ClPl =l\/\PlP2P3 



A \U , ^o fP2P3 P1P2 PlP3\ .^^^. 

V PiP2P3Hi + —\ . 2.22 

2 V Pi P3 P2 J 



The mean Hubble rate H of the mean scale factor a = (010203) ' is given by 



^ = ^^ = > + ^^ + ^=>- P'^^' 



and the Friedmann equation is 

i^^ = ^P+^-^-|ni>), (2.24) 

where the energy density of the scalar field is p = Pt/'^\PiP2P3\, the shear term is given by 

a' = i[(iJi - H,y + {H2 - H,f + {H, - H,)\ (2.25) 

and the potential is 



V{p) = -^— 
P1P2P3 



2 , 2 , 2^ (P2P3\^ fP3PlY fPlP2^^ 



2{pi+p',+pi) 



Pi J \ P2 J \ P3 



(2.26) 



Clearly, these dynamics are quite complex already at the classical level and, as mentioned 
in the introduction, become chaotic as a singularity is approached. The one exception is the 
case when the matter field is a massless scalar field which is precisely what is considered 
here. In this case, as the singularity is approached, the Friedmann equation is asymptotically 
velocity term dominated (AVTD) which means that the potential can be safely neglected [18]. 
Thus, as the singularity is approached, the dynamics are the same as those of the Bianchi I 
space-time with a massless scalar field. This behaviour will be important for the study of the 
effective equations later. However, the quantum Hamiltonian constraint operator derived in 
the following section will hold everywhere and it will be relatively straightforward to extend 
it for other types of matter fields which classically allow the full Mixmaster dynamics. 



Finally, before moving on to the quantum theory, let us consider the parity transformation 
Ilfc which flips the kth physical triad vector e^. (Keep in mind that this transformation does 
not act on any of the fiducial quantities which carry the label o.) These correspond to 
residual discrete gauge transformations. Under this map, we have: Qab -^ Qab, ^abc -^ ^abc 
but eijk —7- —eijk, e — > —e. The canonical variables Ci,pi transform as proper internal vectors 
and co-vectors. For example, 

ni(ci,C2,C3) -> (-Ci,C2,C3) and Ili{pi,p2,P3) ^ {-pi,P2,P3) ■ (2.27) 

Consequently, both the symplectic structure and the Hamiltonian constraint are left invari- 
ant under any of the parity maps 11^ . 

The Hamiltonian description given in this section will serve as the starting point for the 
loop quantization in the next section. 



III. QUANTUM THEORY 

This section is divided into three parts. In the first, we discuss the kinematics of the 
model and in the second we introduce the Hamiltonian constraint operator and describe its 
action on physical states. Finally, in the third subsection we show that the dynamics of a 
wave function sharply peaked around an isotropic geometry are well approximated by the 
LQC dynamics of the closed FRW model. 

A. LQC Kinematics 

The kinematics for the LQC of Bianchi IX models is identical to that of the Bianchi II 
models [16], but we will briefly present the kinematics here as well for the sake of complete- 
ness. 

The elementary functions on the classical phase space that have unambiguous analogs in 
the quantum theory are the momenta pi and holonomies /i^^ of the gravitational connection 
Al^ along the integral curves of el of length fj,io with respect to the fiducial metric qab- These 
holonomies are given by 

hi^\ci, C2, C3) = exp (yUCfcTfc) = cos — ^I + 2 sin ^t/,, (3.1) 

where the r^ are — i/2 times the Pauli matrices. This family of holonomies is completely 
determined by the almost periodic functions exp{ifiCk) of the connection. These almost 
periodic functions will be the elementary configuration variables which will be promoted 
unambiguously to operators in the quantum theory. 

It is simplest to use the ^-representation to specify the gravitational sector 'H|[^^ of the 
kinematic Hilbert space. The basis is orthonormal in the sense that 

{P1,P2,P3\P[,P2,P3) = Sp,p[Sp2P'2^P,,p'^ , (3-2) 

where the right side features Kronecker delta symbols rather than Dirac delta distributions. 
Kinematical states consist of countable linear combinations 

1^) = Yl '^(P^^P^^P3)\Pi^P^^P3) (3-3) 

Pl,P2,Pi 



of these basis states for which the norm 

m\' = E i^(pi.P2,P3)p (3.4) 

pi ,P2 ,P3 

is finite. 

Next, recall that on the classical phase space the three reflections IIj : e^ — )■ — e^ are 
large gauge transformations under which physics does not change since both the metric and 
the extrinsic curvature are left invariant. These large gauge transformations have a natural 
induced action, denoted by IIj, on the space of wave functions \E'(pi,p2,P3). For example, 

ni^(pi,P2,P3) = ^(-Pi,P2,P3). (3.5) 

Since Jlf is the identity, for each i the group of these large gauge transformations is simply 
Z2. As in Yang-Mills theory, physical states belong to its irreducible representation. For 
definiteness, as in the isotropic and the Bianchi type I and type II models, we will work 
with the symmetric representation. It then follows that 'H^i^ is spanned by wave functions 
'^{pi,P2,P3) which satisfy 

^(Pl,P2,P3)=^(bl|,b2|,b3|) (3.6) 

and have a finite norm. 

The action of the elementary operators on 'H^j^ is as follows: the momenta act by mul- 
tiplication whereas the almost periodic functions in q shift the ith argument. For example, 

[pi'^]{pi,P2,P3) =Pi^{Pi,P2,P3) and exp(?;UCi)^ {pi,p2,P3) = ^{pi-87r'jGhiJ,,p2,P3) ■ 

__ __ (3.7) 

The expressions for p2,exp{ifj,C2),P3 and exp^ific^) are analogous. Finally, we must define 
the operator i since e features in the expression of the Hamiltonian constraint. Following 
[16], we define 

, f \Pl,P2,P3) iipiP2P3>0, , , 

£bi,P2,P3) ■■={ . . .^ ^ (3.8) 

[ -bi,P2,P3) if P1P2P3 < 0. 

Finally, the full kinematical Hilbert space "Hkin will be the tensor product "Hkin = ^kin^ ® 
^kin**' where "H^^** = L^(R, dT) is the matter kinematical Hilbert space for the homoge- 
neous scalar field. On "HJJJq**, T will act by multiplication and Pt '■= —thdrp will act by 
differentiation. 



B. The Quantum Hamiltonian Constraint 

To define the quantum Hamiltonian constraint, we must express the Hamiltonian con- 
straint in terms of almost periodic functions of the connection which can be directly pro- 
moted to operators. For isotropic and/or spatially flat space-times, this can be done by 
expressing the fleld strength F^f,^ in terms of holonomies and this is what is done for the 
fii approach in LQC in [8, 10, 11, 15]. However, this is not possible for space-times which 
are both anisotropic and spatially curved such as the Bianchi type II and type IX models. 
In this case we need to extend the strategy: the connection itself — rather than the fleld 



strength — has to be expressed in terms of holonomies. This task was carried out in [16]. 
The connection operator is given by 



where 



/ii 



' bi|A4, 

\P2P3\ 



Ck 



1^2 



sin(/ifcCfc) 
fik 



'b2|A4i 



bipsl 



1^3 



' b3|A4, 

\PlP2\ 



(3.9) 



(3.10) 



and A£|[ = Ay^ivj ip^ is the 'area gap'. Note that the choice for this operator is motivated 
by LQG: it is obtained in [16] by expressing the connection in terms of holonomies, a 
procedure commonly used in LQG, and then ensuring that this approach is equivalent to 
what is done for simpler cosmological models. Although the precise value of the area gap 
may change as the relation between LQG and LQG is better understood, the form of jli in 
terms of the Pi is necessary in order to obtain the correct infrared, low curvature behaviour. 
Using the connection operator, it is possible to promote the classical Hamiltonian con- 
straint in Eq. (2.16) to an operator. Ignoring factor ordering ambiguities and inverse triad 
operators for the moment, Ch is given by 



C 



H 



87rG72A4i 



PiP2\P3\ sin/iiCisin/i2C2 + \p1\p2P3 sm iJ,2C2 sin ^303 



+ Pi \P2 \P3 sin fiscs sin /xiCi 



^qC 



inG-fWAt 



PI 



P1P2 



+ P2P3\ 



P1P2 

P3 



Ipml 
\pl\ 

2 



sm/iici +P3Pi^ 



\P3Pl\ 



P2P3 
Pi 



P2 J 



\P2 
2" 



■Sm/i2C2 







sm /igca 



327rG'7 



2{pI+pI+pI) 



(3.11) 



where for simplicity of notation here and in what follows we have dropped the hats on the 
Pi and sin /ijCj operators. 

To obtain the action of the sin/ijCj operators (or, equivalently, the exp{ifiiCi) operators) 
we will use the same strategy as in [15]. As shown there, it is simplest to introduce the 
dimensionless variables 

sgn(p.)/N (3^^2) 



' (47r7yA4^)i/3 • 

Then the kets |Ai,A2,A3) constitute an orthonormal basis in which the operators pk are 
diagonal 

Pfc|Ai,A2,A3) = [sgn(A,)(47r7v^4i)'^'A^] |Ai,A2,A3), (3.13) 

and quantum states are represented by functions ^(Ai, A2, A3). Then the operator e*^^'^^ acts 
by shifting the wavefunction, 

[e^^i^i M/] (Ai, A2, A3) = ^(Ai - ^^, A2, A3) 



^ 



IA2A3I 
,t;-2sgn(A2A3) 



Ai, A2, A3 



(3.14) 



10 



where we have introduced the variable v = 2A1A2A3 which is proportional to the volume V 
of the space-time: 

V^(Ai,A2,A3) = [27r7yA|t;|4j^(Ai,A2,A3). (3.15) 

The action of the operators e*^^*^^ and e*'^*'^^ is analogous. 

We are now ready to write the Hamiltonian constraint explicitly in the Aj-representation, 
again ignoring factor-ordering issues for the time being: 



Ch = Ci + 62 + 4 + C4 + Ipt, 



(3.16) 



where 



Ci 



-IwMliV^ 



sgn(Ai A2) sin /iiCi sin /i2C2 + sgn(A2 A3) sin /X2C2 sin n^c^ 

+ sgn(A3Ai) sin /igcg sin /iici ; 

1 1 

(AiAs)^— ;==sin;U3C3 + (A2A3)^— ;== sin^uid 

vbsl vbil 

1 



4 = -2nVAh4loe 



C, 



+ (AsAi)^— ==sm/i2C2 

«o(l + 7 ) Ai + A2 + A3 



47 



ll 



C4 = i(167rVA)^/='7rAn4/^(i + ^2) Ux^X2t\ + (AsAs)'^ + (A3Ai)^ 

I Pi Pt Pi 



(3.17) 

(3.18) 

(3.19) 
(3.20) 



It will be straightforward to deal with Ch since the terms in Ci are the exact terms that 
appear in the Bianchi I model and have already been studied in [15] while the terms in C2 
and C4 are of the same form as some of the terms in the Bianchi II model [16] . Finally, the 
only new terms — those in C3 — act by multiplication and will not cause any difficulty. 

All of the terms will be factor-ordered in a symmetric manner. For example, the first 
term in Ci will be factor-ordered as 



"Y^ttMIiv^ 



(sin/iiCiSgnAi + sgnAi sin/iiCi)|T;|(sin/i2C2SgnA2 + sgnA2 sin /i2C2) 



+ (sin/i2C2SgnA2 + sgnA2sin/i2C2)|t;|(sin/iiCisgnAi + sgnAi siuyUiCi) 



vl (3.21) 



while the first term in Co will be 



7rVAM^i£o(AiA2) 



bsl^/^ 



e sin /i3C3 + sin ^^c^e 



bsl^/^' 



(3.22) 



Since all of the components in each term in C3 and C4 commute, there are no factor-ordering 
choices to be made for these terms. 

The factor ordering given in Eq. (3.21) was first introduced in [17] for the study of the 
Gowdy model. It is a particularly nice choice as it causes the octants to decouple from each 
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other, one can then focus on the dynamics of a single octant and then derive the behaviour 
of the other octants via the parity properties of the wave function. 

The only operators that remain to be defined are the inverse volume operators. Using a 
variation on the Thiemann inverse triad identities [29], one obtains the operator [16] 

|prP^|Ai,A2,A3) = '^'^''^^'^-^(^^ (v/|t; + sgn(A2A3)| - Vl^ - sgn(A2A3)|) |Ai,A2,A3). 

(3.23) 
This operator is diagonal in the eigenbasis |Ai, A2, A3) and, on eigenkets with large volume, 
the eigenvalue is indeed well approximated by |pi|~^''^, whence on semi-classical states it 
behaves as the inverse of |p|^/^, just as one would hope. Nonetheless, there are interesting 
nontrivialities in the Planck regime, the most important one being that the inverse triad 
operator annihilates states |Ai, A2, A3) where v = 2A1A2A3 = 0. 

Finally, the other inverse triad operator which is necessary for the study of Bianchi IX 
models can be defined by 

P7' ■■= (ipf^')' ■ (3.24) 



Note that both of these operators were already introduced for the study of the Bianchi II 
model in [16]. 

As in the Bianchi I model, the action simplifies if we replace (Aj, Aj, A^) by (Aj, Aj,w)^. 
Because of the high symmetry of the Bianchi IX model, it does not matter which of the 
Aj is replaced; we will choose to replace A3 by v here. This change of variables would be 
nontrivial if, as in the Wheeler-DeWitt theory, we had used the Lesbegue measure in the 
gravitational sector. However, it is quite tame here because the norms are defined using a 
discrete measure. The inner product on "Hf^^^ is now given by 

(^i|^2)kin= Yl '^i(Ai,A2,t;)^2(Ai,A2,t;) (3.25) 

Ai,A2,D 

and states are symmetric under the action of 11^. In the Appendix of [16], it is shown that 
under the action of the IIj, the operators sin/ijCj have the same transformation properties 
as Cj under the reflections IIj in the classical theory. As a consequence, Ch is also reflection 
symmetric^. Therefore, its action is well defined on "Hf^^^: Ch is a densely defined, symmetric 
operator on this Hilbert space. In the isotropic and Bianchi I cases, its analog has been shown 
to be essentially self-adjoint [30, 31]. In what follows we will assume that (3.16) is essentially 
self-adjoint on "Hf^^^ and work with its self-adjoint extension. 

We can now study the action of Ch on a wavefunction. For a complete derivation of the 
action of each term in the constraint, see [15, 16]. 

It is straightforward to write down the full Hamiltonian constraint on "Hl^^^: 

-h''d^'^{X2,X3,v;T) = e^{X2,X3,v;T), where 6 = -2Cgrav. (3.26) 



^ This cannot be done for states where A1A2A3 — but since these states decouple under the action oiCn, 
we can restrict our attention solely to states where A1A2A3 7^ 0. 

^ Note that although Xlielli — —e (recall that classically e -^ — e under a parity transformation) only when 
11 7^ 0, in the v = case the wavefunction is annihilated by the gravitational part of the Hamiltonian 
constraint Cgrav and therefore HiCgravHi losing) = = Cgrav losing) where losing) is a state that only has 
support on V = 0. It is then straightforward to show that tliCniiil'^) = Ci/|^) for all wavefunctions. 
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As in the isotropic case [32], one can obtain the physical Hilbert space T-iphy by a group 
averaging procedure and the final result is completely analogous. Elements of 'Hphy consist 
of 'positive frequency' solutions to (3.26), i.e., solutions to 



-thdT^{\i,\2,v;T) = ./\e\^i\u\2,v;T), 
which are symmetric under the three refiection maps IIj: 

^(Ai,A2,t;;T) = ^(|Ai|,|A2|,|t;|;T). 
The scalar product is simply given by 

(^i|^2)phys = (^i(Ai,A2,t;;To)|^2(Ai,A2,t;;To))kin 

Ai,A2,?J 



(3.27) 



(3.28) 



(3.29) 



where To is any "instant" of internal time T. 

Since elements of "H^^^^ are invariant under the three parity maps 11^ and the Hamiltonian 
constraint satisfies flfcCgravn/fc = Cg^av, knowledge of the restriction of the image Cg^av^ of ^ 
to the positive octant suffices to determine Cgrav^ completely. Therefore, in the remainder 
of this section we will restrict the argument of C//^ to the positive octant. The full action 
is simply given by 

(Cg,av^)(Al,A2,t;) = (Cgrav^)(|Al|,|A2|,|t;|). (3.30) 



Since all states with v = are annihilated by Cgrav, their evolution is trivial: 

4^(Ai,A2,t; = 0;T) = 0. 



(3.31) 



Such states correspond to classical geometries which are singular and therefore we will call 
these states 'singular', even though they are well defined in the quantum theory. Non- 
singular states on the other hand are physically much more interesting. On them, the 
explicit form of the full constraint is given by: 



d^^{\i,\2,v;T)=nG 



("(t; + 2)v/^7T4 ^+(Ai, A2, t;; T) - (t; + 2)v/^^+(Ai, A2, t;; T) 



- 9,^2iv - 2)v^^o-(Ai, \2,v; T) + O.^^iv - 2) ^1^; - 4| ^4 (Ai, \2,v; T) 
'^'^''^-, (v/^^TT - v^^^^) ($+ - B.-2^-\ (Ai, A2, V, T) 



(47r7VA) 
+ ^.^^°^^^.T!^ f (v/^^TT - v^^^^)^ (AiA2)^ + (A2A3)« + (AsAO^ 



(47r7\/A)2/3 
-l(At + A^ + A^)')M/(Ai,A2,t;;T) 



where dx is the step function 



1 if X > 0, 
if a; < 0. 



(3.32) 



(3.33) 
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Note that the step function kills any terms that would allow the positive octant to interact 
with any of the other ones, this is a direct consequence of the factor ordering choices made 
earlier. 

The ^^^4 are defined as follows: 

/q; 4- A 7; + 2 \ /?; ± 4 

vI/f(Ai,A2,t;;T) = v[/ -Ai, \2,v ±4-t) +^ (^ ■ \„\„v ±4-T 

\v ±2 V J \v ±2 

+ vtf!i±l.A„^.A„.±4;r)+«('^:±l.A,.A„.±4;r 

\ V V ±2 J \ V 

+ ^ (Xi, ^ ■\2,v±4;t]+^ (x,, ^ ■\2,v± 4;t] , (3.34) 

and 

v^o^(Ai, A2, v;T)=^i! f ^^ ■ Ai, ^ ■ A2, v;t)+^ f ^^ ■ Ai, A2, v; T 
\ V V ±2 J \ V 

[ V V ±2 \ f V 

+ ^ —— ■ Ai, A2, t;; T + ^ —— ■ Ai, A2, v- T 

\v ± 2 V J \v ± 2 

+ ^ (ai, ^ ■ \2,v;t\ + ^ (Xu ^ ■ \2,v;t\ , (3.35) 

while ^"^ are given by 

$±(Ai,A2,t;;T) =(^1^; ± 2 + 1| - ^1^^ ± 2 - 1|) x (A2A3)'^(^ ■ Ai, A2,t; ± 2; T 

+ (A3Ai)^^(Ai, ^ ■ A2,^ ± 2; t) + (AiA2)^^(Ai, A2, t; ± 2; T 

(3.36) 



As expected, the quantum dynamics of the Bianchi IX model reduces to that of the 
Bianchi I model discussed in [15] in the limit £0 — )■ in Eq. (3.32). 

Eq. (3.32) also immediately implies that the steps in v are uniform: the argument of the 
wave function only involves f — 4, f — 2, f , f + 2 and v + 4. Thus, there is a superselection 
in V. For each e G [0,2), we can introduce a lattice C^ of points v = 2n + e. Then the 
quantum evolution — and the action of the Dirac observables pT and V\t commonly used in 
LQC — preserves the subspaces "Hp^y consisting of states with support in v on C^. The most 
interesting lattice is the one corresponding to e = since it includes the classically singular 
points V = 0. 

The form of the action of the Hamailtonian constraint operator also shows that the 
classical singularity is resolved. Using the scalar field T as time, we find that if one starts 
with a wavefunction which only has support on singular states, that wavefunction does not 
evolve in T and therefore will always only have support on singular states. 

On the other hand, a state which does not have any support on the singular subspace 
will never have support on it. Restricting our argument to the positive octant for the sake 
of simplicity (it can easily be generalized to the other octants), it is easy to see that to go 
from Ai, A2, f > to f = 0, one must either have v = 2 and then $~ will give a term with 
w = or have v = 4 and then \1/J will give a term with v = 0. However, the prefactors in 
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front of $~ vanish for f = 2 just as the prefactors in front of "^^ vanish for v = 4. Because 
of this, it is impossible for a wavefunction with no support on singular states to ever gain 
support on a singular state. 

This shows that singular states decouple from nonsingular states under the relational T 
dynamics given by Eqs. (3.31) and (3.32). In other words, if one starts with a nonsingular 
state at some 'time' To, it will remain nonsingular throughout its evolution. It is in this 
(rather strong) sense that the singularity is resolved. 



IV. EFFECTIVE EQUATIONS 

In the isotropic models, effective equations have been introduced via two different ap- 
proaches — the embedding [33, 34] and the moment expansion [35] methods — in order to 
study the first order quantum-corrected equations of motion. In the isotropic case the effec- 
tive equations following from the embedding approach provide an excellent approximation 
to the full quantum evolution of states which are Gaussians at late times, even in the A 7^ 
as well as k=±l cases where the models are not exactly soluble. However, the truncation 
method is more systematic and also more general in the sense that it is applicable to a wide 
variety of states. Nonetheless, in this section we will use the first method (although we 
will ignore the effect of fluctutations in this work) in order to gain qualitative insights into 
modifications of the equations of motion due to quantum geometry effects. 

To obtain the effective equations we can restrict our attention to the positive octant of the 
classical phase space (where e = 1) without loss of generality. Then the quantum corrected 
Hamiltonian constraint is given by the classical analogue of (3.11): 



<-) ' "-'grav 



7Cff 



0, 



(4.1) 



where^ 



c: 



eff 
grav 



P1P2P3 

87rG'72A4i 



sin /iiCi sin jj,2C2 + sin JI2C2 sin fx^c^ + sin fi^c^ sin /iiCi 



{P1P2) 



3/2 



^' . 2 



sm /X3C3 



{P2P3) 



3/2 



sm/iiCi 



(PsPi) 



3/2 



327rG'72 



2{pI+pI+pI) 



P1P2 

P3 



P2P3 
Pi 



■ sm H2C2 



P3P1 



P2 



(4.2) 



Using the expressions (3.10) of /i^, it is easy to verify that far away from the classical 
singularity — more precisely in the regime in which the Hubbies rates Hi are well below 
the Planck scale — the effective Hamiltonian constraint (4.1) is well-approximated by the 
classical one given in Eq. (2.16). 



^ Recall that every ig which appears in the constraint is divided by Vg which has been set to 2. As Iq/to is 
dimensionless, we must ignore £0 when counting units. 
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The effective dynamics are obtained by taking Poisson brackets with the effective Hamil- 
tonian constraint. This gives 



Pi = 7"^ — (sin/i2C2 + sin /igCg) 



0P2P3 cos/iiCi; 



(4.3) 



Cl 



P2P3 

A4i 



sin /iiCi sin /i2C2 + sin fiiCi sin fi^c^ + sin /i2C2 sin /isCs 



H ;^ COS/iiCi(Sin/X2C2 + Sm/isCsj -— COS/i2C2(,Sin/iiCi +Sin/i3C3j 



2 
/U3C3 



COS;U3C3(Sin/iiCi +Sinfl2C2, 



3 

2^ 



P1P2 . _ , P1P3 . _ 
sm /i3C3 H sm /i2C2 



P3 



P2 



P2P3 . _ 

- — sm^uici 
3pi 



+ 7T CiCOSyUiCi - -P2C3C0S;U3C3 - -P3C2 COS ;U2C2 

2 Pi 2 2 



+ ^(1 + 7^) ( 4pi - 2pi f 4 + 4 V 2#) 
4 V VPs P2/ Pi / 



(4.4) 



The equations for P2,P3, C2 and C3 are the same modulo the appropriate permutations. Note 
that it is easy to extend this for other matter fields and also to the vacuum case simply by 
appropriately modifying the matter part of the effective Hamiltonian constraint. 

In the embedding approach these effective equations provide quantum geometry correc- 
tions to the classical equations of motion Eqs. (2.18) and (2.19) due to the area gap. However, 
careful numerical work comparing the full quantum dynamics to the effective dynamics is 
necessary to determine whether the effective equations are accurate beyond first order in h. 

Now, it is well known that classical Bianchi IX space-times with a massless scalar field as a 
matter source behave in an asymptotically velocity term dominated (AVTD) manner^, that 
is to say that the potential term is negligible (see [18] and references therein). For certain 
regions of phase space, this will occur before quantum gravity effects become important and 
we will assume that in this case only quantum gravity corrections to the velocity terms are 
relevant. 

It then follows that this behaviour is identical to that of the Bianchi I model and therefore 
the effective Friedmann equation for the Planck regime to first order in h is given by [36] 



H' 



871G 



1- 



P 



S2 



■^2\2 



SV 

2pc 327tGp, 



+ o{e 



Pl); 



(4.5) 



where pc = S/Sn-f'^AGil^ ^ 0.41ppi (recall that A = 4^3717 and 7 
hole entropy calculations [37]). The expression for S^ is given by 



0.2375 due to black 



(piCi - P2C2)^ + (P2C2 - PsCsf + (P3C3 - PlCi) 



3^2p3 

and one can show that p^S^ is a constant in the AVTD limit [36]. 



(4.6) 



^ This is true so long as the constant of motion p^ is large enough so that the three scale factors are all 
decreasing as the singularity is approached. 
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It is clear that there is a bounce {H^ = 0) when the matter energy density reaches 

(4.7) 



1 

Pbouncc 



2 



3S2 // 3S2 \ 




Pc - TTT^ + \ \Pc- TTT^ P' 



IQttG V 167rG/ V levrG 






at which point the energy density and curvature will both decrease and leave the Planck 
regime and the classical dynamics will once again become a good approximation. It follows 
that the matter energy density is always bounded above by the critical energy density 
Pc = 0.41ppi. This is only an upper bound as the matter density at the bounce depends 
quite strongly on S^ which is a measure of the strength of the gravitational waves: the 
stronger the gravitational waves are, the lower Pbounce will be. 

The scenario described above relies on the AVTD behaviour of the Bianchi IX cosmology 
with a massless scalar field occuring before quantum gravity effects become important. In 
this case, the true Friedmann equation can then be well approximated by Eq. (2.24) in the 
classical regime and by Eq. (4.5) in the AVTD limit. However, this scenario will not be 
valid for all regions of phase space, in particular where the scalar field momentum pt is 
small enough for the chaotic Mixmaster behaviour to appear. 

It has been suggested that, by bounding the strength of the potential terms due to inverse 
triad effects, quantum gravity effects could play an important role in Bianchi IX dynamics 
and that the chaotic Mixmaster behaviour would be avoided as a result of this for all types 
of matter fields [24] . In the effective equations presented above, we have ignored the effect 
of inverse volume corrections (which for the inverse volume operator used in this paper are 
only important for f < 4) and have only considered the effect of holonomy corrections. If 
the chaotic behaviour is to be generically avoided in this effective theory, it will be because 
the repulsive quantum gravity effects will ensure that the Bianchi IX space-time will not 
remain in high curvature regions for long enough for there to occur a sufficient number of 
Mixmaster bounces for chaos to appear. 

For now, this remains a conjecture and one would have to study the Bianchi IX effective 
equations of motion more carefully, using both analytic and numerical methods, in order to 
determine whether the bounce is generic and also to see if chaotic behaviour is avoided or 
not in the effective theory for small p^. 

V. DISCUSSION 

In this paper we have studied the improved LQC dynamics for Bianchi IX cosmologies 
where the matter content is a massless scalar field which is used as a relational time param- 
eter. We have shown that the singularities in the classical theory are resolved by quantum 
gravity effects in the usual manner in LQC as the singular states decouple from the regular 
ones under the relational dynamics given by the Hamiltonian constraint operator. 

It is important to point out that all of the tools necessary for the task of deriving the LQC 
dynamics for Bianchi IX models were already available. First, the form of fli was introduced 
in the study of Bianchi I models [15], as were the variables Aj which greatly simplify the form 
of the action of the Hamiltonian constraint operator. The other two necessary ingredients 
to the results for this work are the connection operator and the inverse triad operators, both 
of which were introduced for the study of Bianchi II models in [16]. In addition, even the 
fact or- ordering choices necessary in the Hamiltonian constraint operator had been made in 
[15-17]. Because of this, it is reasonable to expect that no additional machinery should 
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be necessary in order to complete the study of the loop quantum cosmology of the other 
Bianchi models of type A. 

Finally, in addition to obtaining a well-defined LQC Hamiltonian constraint operator for 
Bianchi IX space-times and studying some of its properties, we also derived some effective 
equations which provide modifications to the classical equations of motion due to the area 
gap which is a manifestation of quantum geometry in LQG. Although all of the results 
presented in this paper were derived for the particular case of a massless scalar field as the 
matter field, it will be easy to extend the results presented here for other types of matter 
fields (as well as the vacuum case) for both the quantum and effective theories. 

Of course, it is not enough to know the form of the equations of motion given in Eq. (3.32) 
in order to understand the full dynamics of the loop quantum cosmology of Bianchi IX 
models. Numerical studies will be particularly useful and help us understand how the 
quantum state of a Bianchi IX cosmology evolves with time. Most interesting would be a 
study of states which are sharply peaked around a semi-classical state at late times and to 
then evolve them back in time to see what happens as the curvature increases. Based upon 
previous experience with isotropic models, one might expect to see one or several bounces 
as the curvature reaches the Planck scale but careful numerical studies are needed to check 
this. 

If the BKL conjecture is correct a good understanding of the quantum dynamics of Bianchi 
IX cosmologies will lead to a better understanding of the behaviour of generic space-times as 
their curvature reaches the Planck scale. If Bianchi IX models are sufficiently rich in order 
to understand the approach to such regions, it would appear that no singularities would 
form since an initially nonsingular Bianchi IX wave function must remain nonsingular as 
shown in Sec. III. It is therefore possible that a careful study of the BKL conjecture at the 
level of the quantum dynamics could provide a no-singularity theorem, a first step in this 
direction is provided by [22]. 

A simpler avenue to study quantum gravity effects in Bianchi IX models would be to 
study the effective dynamics presented in Sec. IV. In isotropic models it turns out that the 
effective dynamics are surprisingly accurate even in the deep Planck regime: the effective 
equations accurately predict the quantum trajectory throughout the quantum bounce for 
sharply peaked wave functions. Because of this, it would be interesting to study the dynamics 
given by the effective equations for Bianchi IX space-times. However, it is essential to see 
where the effective equations break down, if they do at all. This can be done by including 
higher order corrections to the effective equations via the moment expansion method and 
also by comparing the predictions of the effective equations to full numerical solutions of 
the Hamiltonian constraint operator. 

An analysis of the effective equations of motion in the case where the asymptotically 
velocity term dominated behaviour begins before quantum gravity effects become important 
shows that there is a bounce when the curvature reaches the Planck scale and that the 
matter energy density is bounded above by the critical energy density pc ~ 0.41ppi. This 
result relies on the AVTD behaviour and is not generic. Therefore, one must also examine 
other areas in the phase space in order to fully understand the predictions of the effective 
theory, particularly near Planck scales. As the Mixmaster behaviour appears for small p\ 
during the approach to the singularity in the classical theory, the effective equations can 
provide a better understanding of how quantum gravity effects may modify the Mixmaster 
behaviour as well. In particular it is possible that, as for simpler isotropic models and 
in AVTD case here, these quantum gravity effects will be repulsive and cause a quantum 



18 

bounce. This would limit the amount of time that the Mixmaster behaviour occurs and the 
chaos which arises in the classical theory might be avoided due to the short time span of 
the Mixmaster dynamics. However, this remains a conjecture and much more work, both 
analytic and numerical, is needed in order to resolve this question. 

Finally, it has been pointed out that the dipole cosmology model can be used in order to 
study the Bianchi IX model [25] . Although that paper studies the Euclidean theory, it would 
nonetheless be interesting to compare the model presented in [25] with the one developed 
in this paper. In particular, [25] suggests two possible approaches in order to obtain the 
Hamiltonian constraint operator for their model. Comparing the quantum dynamics result- 
ing from these two possibilities to those derived in this paper could help determine which of 
the two approaches is the correct one and hence give some insight into the dipole cosmology 
models and also spin foam models in general. It is also possible to further probe the relation 
between the canonical and the covariant approaches to LQG via LQC by extending the 
Feynman path integral construction given in [38] for the fiat FRW model to the Bianchi IX 
model; this extension would be nontrivial due to the additional degrees of freedom present 
in Bianchi IX space-times, but it could also improve our understanding of the connection 
between the canonical and covariant approaches to LQG as well as the relation between full 
LQG and the symmetry-reduced models of LQC. 
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